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ABSTRACT: Search for a proper and realistic equation of state (EOS) for strongly interacting 
matter used in the study of the QCD phase diagram still appears as a challenging problem. 
Recently, we constructed a hybrid model description for the quark gluon plasma (QGP) as 
well as hadron gas (HG) phases where we used an excluded volume model for HG and a 
thermodynamically consistent quasiparticle model for the QGP phase. The hybrid model 
suitably describes the recent lattice results of various thermodynamical as well as transport 
properties of the QCD matter at zero baryon chemical potential (/is). In this paper, we 
extend our investigations further in obtaining the properties of QCD matter at finite value 
of jiB and compare our results with the most recent results of lattice QCD calculation. 
Finally we demonstrate the existence of two different limiting energy regimes and propose 
that the connection point of these two limiting regimes would foretell the existence of critical 
point (CP) of the deconfining phase transition from HG to QGP in our hybrid model. 
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1 Introduction 

Quantum chromodynamics (QCD) predicts that at sufficiently high temperatures (T) and/or 
chemical potentials strongly interacting matter goes through a phase transition from 

colour insulating hadron gas (HG) phase to colour conducting quark gluon plasma (QGP) 
phase [1, 2]. Heavy ion collisions provide a unique opportunity to study this QCD phase 
transition from HG to QGP. A remarkable outcome of these experiments is the observation 
of a hot, strongly interacting matter which behaves much like a nearly perfect relativistic 
fluid [3]. Our endeavour in such studies is the search of a suitable equation of state (EOS) 
for the description of both phases of strongly interacting matter. Significant success has 
been gained in lattice calculations using QCD thermodynamics to provide a valid EOS for 
QCD matter at zero baryon chemical potential [4-9]. However, the lattice methods still 
lack a reliability to describe the properties of matter possessing a finite density of baryons. 
Therefore, finding an EOS for QCD matter valid at zero as well as non-zero chemical po- 
tential is still a challenging problem. Consequently, in the absence of a proper and realistic 
EOS for QCD matter, precise mapping of the entire QCD phase boundary and finding the 
location of a hypothesized QCD critical point (CP) appears as a tough challenge before the 
experimentalists and theorists [10-12]. 

Recently, we constructed a hybrid model for HG and QGP phases where we made 
use of a new excluded volume model for HG and a thermodynamically consistent quasi- 
particle model for the QGP phase [13, 14]. We have compared various thermodynamical 
and transport properties at fj,s = as obtained in our hybrid model description with the 
recent lattice results and found an excellent agreement between them [13, 14]. We have 
also constructed a deconfining first order phase boundary between HG and QGP and found 
an end point where the boundary terminates [14, 15]. Further, we have shown that at 
this end point, the nature of phase transition changes as As/T 3 = (sqgp — shg)/T 3 and 
Ac 2 , = (c 2 )qgp — {c 2 s )hg become zero where s is the entropy density and c 2 is the square 
of speed of sound [13]. We have also found a cusp like behaviour in the variation of rj/s 
(where rj is the shear viscosity) with respect to temperature at the end point obtained in 
our hybrid model [13] which was, therefore, suggested as a hint for the existence of QCD 
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CP by some authors [16, 17] earlier. All these results suggest that this end point could be 
taken as a CP of deconfinement phase transition from HG to QGP. 

In this paper, we want to explore further the suitability of the hybrid model in predicting 
the properties of QCD matter having a finite value of As we have pointed out that a hot, 
strongly coupled and nearly perfect fluid was observed at RHIC, thus ideal fluid description 
can indeed be appropriate to explain the expansion history of the matter created during the 
heavy ion collisions [18, 19] at least at RHIC and LHC. Ideal relativistic hydrodynamics 
is entropy conserving i.e., the ratio of the entropy and the particle number is expected 
to remain constant during expansion of the medium in heavy ion collision. The EOS of 
strongly interacting matter should indicate such hydrodynamic description of the medium. 
Very recently Wuppartel-Budapest (WB) collaboration gave the isentropic trajectories on 
the T — Hb plane and calculated various thermodynamical as well as transport properties 
along these trajectories [20]. This lattice calculation involves Taylor's expansion method 
to circumvent the fermion sign problem which always exists at finite [Ib- In the same 
paper [20], authors have also shown the results for various thermodynamic quantities at 
Hl = 400 MeV and fig = 400 MeV (where //b/3 = fj, u = y,^ = y s and = \x u = 

while strangeness chemical potential fj, s has been calculated by imposing the net strangeness 
neutrality condition i.e. n s = 0). Therefore, it will be worthwhile to compare the results 
obtained by the collaboration with the corresponding results of our hybrid model and thus 
its suitability for the description of QCD matter can be tested. 

We must emphasize the new results obtained by us in order to demonstrate the impor- 
tance of the paper. We first draw the isentropic trajectories at some fixed S/N (which is 
equivalent to fixed S/Nl trajectories of Ref [20] since in our model we include the condition 
of net strangeness neutrality) on T — [Xb plane. We plot the isentropic trajectories having 
S/N = 300, 100, 45 and 30. Here S/N = 300, 45, 30 correspond to typical ratios 
at RHIC, SPS and AGS, respectively [21]. We also show the deconfining phase boundary 
between HG and QGP and the location of CP on this phase boundary as obtained earlier 
in our hybrid model [14, 15]. Further, we plot the chemical freezeout curve as obtained in 
our exclude volume model of HG [22] on this T — hb plane to demonstrate its proximity 
to CP. After it, we calculate various thermodynamic as well as transport quantities like 
normalized pressure density (p/T 4 ), normalized energy density (e/T 4 ), normalized entropy 
density (s/T 3 ), trace anomaly (e — 3p/T 4 ) and square of speed of sound (c^) and plot 
their variations with temperature along an isentropic trajectory having S/N = 300 which 
corresponds to the highest RHIC energy i.e., ^/sjvtv = 200 GeV. We compare our model 
results with the corresponding results obtained in recent lattice calculation [20]. We also 
calculate p/T 4 , e/T 4 and s/T 3 for strongly interacting matter at fJLs = 400 MeV and show 
their variations with respect to temperature in the present hybrid model description which 
is again compared with lattice data. Finally, we calculate a ratio of total number density of 
baryons (ng) to the total number density (tit) and show its value at CP as obtained in our 
model. Based on this ratio we proposed two different limiting regimes of HG as suggested 
earlier [23, 24] and discuss the nature of the phase transition between these two regimes. 
We have also discussed the possibility that the CP results as the meeting point of these two 
limiting regimes and HG thermodynamics shows a clear change at the CP. 
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The rest of the paper is organised as follows : in section 2 and 3, we will provide a 
brief description of our HG model and thermodynamically consistent quasiparticle model 
of QGP, respectively. Detailed description of these models can be found in our previous 
work [13, 14, 25]. Section 4 will demonstrate all the results obtained by us in our hybrid 
model and their comparison with the lattice results. At the end of section 4, we will also 
provide some insights regarding the existence of CP on the deconfining first order phase 
boundary between HG and QGP phases by calculating a ratio n^/ny at CP. 



2 Formulation of excluded volume model for HG 

In the excluded volume models, repulsive interaction between two baryons has been in- 
cluded by giving the baryons a hard-core geometrical size and consequently reduces the 
hadronic degrees of freedom at large T and/or Consequently, the hadronic pressure 
is reduced and one can get a deconfinement phase transition from HG to QGP by using 
Gibbs' construction. We incorporate excluded-volume correction arising due to baryonic 
size only in our excluded volume model. We assume that mesons can overlap and fuse 
into one another and hence do not possess any hard-core repulsion. The grand canonical 
partition function for the HG, with full quantum statistics and after incorporating excluded 
volume correction can be explicitly written as [13, 22, 25]: 
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where gi is the degeneracy factor of ith species of baryons, Ei is the energy of the particle 
(Ei = y k 2 + m 2 ), V® is the eigenvolume of one baryon of ith species and NjV® is the 
total occupied volume by the baryons and Nj represents total number of baryons of jth 
species. 

Now we can write eq.(2.1) as: 
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where Jj represents the integral: 
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and Aj = exp{^) is the fugacity of the particle, n!p is the number density of jth type of 
baryons after excluded volume correction and can be obtained from eq.(2.2) as: 
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This leads to a transcendental equation as 

nf = (1 - R)^ - h\ 2 — + A 4 2 (l - R)k (2.5) 

where l[ is the partial derivative of Jj with respect to Aj and i? = £^ n t x V^ 1S the fractional 
occupied volume. We can write R in an operator equation as follows [25, 26]: 

R = Ri + ClR (2.6) 

where i?i = ^ with R° = £ n?t? + £ 

is the density of pointlike baryons of 

ith species and the operator f2 has the form : 

Using Neumann iteration method and retaining the series upto Cl 2 term, we get 

R = Rx + ClRi + n 2 Ri (2.8) 

Eq. (2.8) can be solved numerically for R. Finally, we get the total pressure [25, 26] of the 
hadron gas: 

pf G = T(l -R)J2 k\i + E P™ ( 2 - 9 ) 

» 3 

In eq. (2.9), the first term on the right hand side represents the pressure due to all 
types of baryons and the second term gives the total pressure from all mesons in HG. In 
this calculation, we have taken an equal volume V° = for each baryon with a hard- 
core radius r = 0.8 fm. The reason to choose this hard-core radius has been discussed in 
ref. [13]. We have taken all baryons and mesons and their resonances having masses upto 
2GeV/c 2 in our calculation of HG pressure. We have also used the condition of strangeness 
neutrality by putting Y^i Si( n f ~ ^| ) =0, where Si is the strangeness quantum number of 
the ith hadron, and nf (nf) ^ s ^ ne s t ran g e (anti-strange) hadron density, respectively. 

3 Quasiparticle model (QPM) 

The EOS for QGP, as used in this paper has been described in detail in reference [27, 28]. 
The effective mass of the gluon in the quasiparticle model is [13, 14, 27, 28]: 

™l(T) = f 9 2 (T)T 2 (l + ^j, (3.1) 

where N c represents the number of colours. We have taken N c = 3 in our calculation and: 

* = + < 3 - 2 > 
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Here Nf is the number of flavours of quarks and ///is the quark chemical potential belonging 
to the flavour f . We take the effective mass of the quarks as [34] : 



m 



qO 



+ V2m q0 m th + ml, 



(3.3) 



were m q o is the rest mass of the quarks. nith represents the thermal mass of the quarks [13, 
14, 27, 28]: 



N 2 



8N, 



T 2 + 4 



(3.4) 



Taking these values for the effective masses, energy density can be derived from the grand 
canonical partition function in a thermodynamically consistent manner and is given as 
[13, 14]: 
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(3.5) 



with 6i(xil) = (xil) 3 Ki(xil) + 3(xil) 2 K2(xil), where K\ and are the modified Bessel 
functions with x% = ^ and index i runs for gluons, up-down quarks q, and strange quark s. 
Here di are the degeneracies associated with the internal degrees of freedom. Now, by using 
the thermodynamic relation e = — p, pressure of system at \x q = can be obtained as 
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where po is the pressure at a reference temperature To. We have used po=0 at To=100 
MeV in our calculation. We get the pressure for a system at finite [i q as : 



v{T,Hq) = p{T,0) + n q d\i q . 
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where the expression for n q is: 
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with Ii(xil) = (xil) 2 K2(xil). Thus all the thermodynamical quantities can be obtained in 
a consistent way by using this model. We have used the following expression for the T and 
Hb dependence of the coupling constant [13] 
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where At = UbMeV and a 
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Figure 1. Isentropic trajectories in T — /j,b plane. First order deconfining phase boundary is 
shown by solid curve and solid diamond symbol represents the location of CP on this deconfining 
phase boundary between HG and QGP. Long dashed curve with solid circular symbols presents the 
chemical freeze-out curve obtained in our excluded volume model of HG. 



4 Results and Discussion 

Figure 1 shows the trajectories for a given S/N (= S/Nl of Ref. [20]) in T — fis plane 
where S is the total entropy, N is the net particle number and iVf, is the net light particle 
number [20]. Here N is equal to iV^ since we include net strangeness neutrality condition 
to calculate fig- We also plot the first order deconfining phase boundary between HG and 
QGP phase as obtained in reference [14, 15] . The end point of this phase boundary is 
shown by solid diamond symbol. We have shown earlier that the phase transition changes 
its nature from first order to higher order [13] at this end point and thus it represents a CP 
of deconfmement phase transition between HG and QGP phases. Furthermore, we draw 
the chemical freezeout curve as obtained in our excluded volume model for HG by fitting 
the experimental data of particle ratios [22]. Our CP corresponds to a value of S/N = 60 
and it corresponds to a freezeout point at T = 158 MeV and = 170 MeV which further 
gives y^s/vTV ~ 21 GeV. From figure 1, it is clear that moving towards higher S/N values 
means moving the trajectories from high region to low region. Here S/N = 300 
thus corresponds to highest RHIC energy i.e., ^/snn = 200 GeV. 

Figure 2 demonstrates the variation of normalized pressure density (p/T 4 ) with respect 
to temperature for a trajectory at a constant S/N (= 300) as obtained from our hybrid 
model. This trajectory corresponds to ^/snn = 200 GeV as mentioned earlier. We also 
compare our results with the corresponding results obtained in a recent lattice calcula- 
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Figure 2. Variation of normalized pressure density with respect to temperature at constant 
S/N = S/N L = 300 corresponding to RHIC energy (i.e. = 200 GeV)[21]. Solid squares 

represent the recent lattice results [20]. 



tion [20] for the trajectory having the same value of S/N. We find that our hybrid model 
provides an excellent fit to the lattice data in the entire temperature range. 

Similarly, figure 3 and figure 4 present the variations of normalized energy density 
(e/T 4 ) and normalized entropy density (s/T 3 ), respectively with respect to temperature. 
We again compare our results with the corresponding results obtained in the lattice calcu- 
lation [20] and find excellent agreement between them. We also observe that in our hybrid 
model result, there is no sign of any discontinuity at the joining point of HG and QGP 
phases in the variations of e/T 4 and s/T 3 with respect to temperatures. It is due to the 
passing of this constant S/N trajectory from HG to QGP phase in crossover region as 
shown in figure 1 in our hybrid model. Thus our model supports a smooth crossover from 
HG to QGP or vice-versa at highest RHIC energy and there is no hint for the occurrence 
of first order phase transition. 

In figure 5, we show the variation of trace anomaly (e — 3p/T 4 ) with respect to temper- 
ature in our hybrid model calculation. This quantity provides a measure of the interactions 
present in the medium. We also plot the lattice data obtained in [20] for a comparison 
with our model results. We find that our hybrid model again appropriately describes the 
lattice data from very low to high temperatures. This study shows that the interaction 
present in the medium produced at highest RHIC energy is quite strong at a temperature 
corresponding to 170 MeV. At higher temperature side, trace anomaly slowly decreases with 
the increase in temperature, however, it remains significant even upto 500 MeV. Therefore, 




Figure 3. Variation of normalized energy density with respect to temperature at constant 
S/N — S/Nl — 300. Solid squares represent the recent lattice results[20]. 
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Figure 4. Variation of normalized entropy density with respect to temperature at constant 
S/N = S/Nl = 300. Solid squares represent the recent lattice results[20]. 
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Figure 5. Variation of trace anomaly factor (e — 3p/T 4 ) with respect to temperature at constant 
S/N — S/Nl — 300. Solid squares represent the recent lattice results[20]. 
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Figure 6. Variation of square of speed of sound (cf) with respect to temperature at constant 
S/N = S/Nl — 300. Solid squares represent the recent lattice results[20]. 
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Figure 7. Variation of normalized pressure density with respect to temperature at /j,b = 400 MeV 
and comparison with lattice results[20]. 




(s/T 4 ) 



LATTICE DATA 



(s/T 4 ) hg ((I b =400 MeV) Our Model 



(s/T 4 ) QGp (n B =400 MeV) Our Model 



50 100 150 200 250 300 350 400 450 500 

T (MeV) 



Figure 8. Variation of normalized energy density with respect to temperature at [is = 400 MeV 
and comparison with lattice results[20]. 
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Figure 9. Variation of normalized entropy density with respect to temperature at /ig — 400 MeV 
and comparison with lattice results[20]. 

inclusion of interactions appears important in deriving the properties of the strongly in- 
teracting matter created at RHIC highest energy. Experimentalists have also found the 
signatures of strongly interacting QGP (sQGP) by measuring the elliptic flow at RHIC [3]. 

Figure 6 presents the variations of square of speed of sound (c^) with respect to temper- 
ature at a constant S/N (= 300) obtained in our model. We find that it suitably matches 
with the Cg obtained in lattice calculation [20]. We observe a shallow minimum in (? s around 
T ~ 160 — 170 MeV which corresponds to a smooth crossover from HG to QGP phase or 
vice- versa as suggested in Ref. [29, 30]. All the results at constant S/N = 300 suggest that 
the medium goes through a smooth crossover from HG to QGP at highest RHIC energy 
i.e., y^STVTV = 200 GeV around T = 170 MeV and the strength of interaction present in 
the medium is still quite sigificant. 

Now, we want to calculate some thermodynamical quantities at a large value of \xb 
(= 400 MeV) and make comparisons with the corresponding results obtained in lattice 
calculation [20]. However, before making this exercise it is very important to note that for 
a large value of /^b, we need other lattice calculations separately in order to support the 
trends given here as the notorious sign problem creeping into lattice calculation and make it 
difficult at net baryon density (or baryon chemical potential). In figure 7, we have presented 
the variations of the normalized pressure density with respect to temperature at /j>b = 400 
MeV. We have compared it with the lattice data points and found a good agreement with 
the results of our hybrid model in the low temperature region i.e., upto 225 MeV. However, 
at higher temperatures our hybrid model results are slightly higher than the lattice points. 
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Table 1. Ratio of tlb/tit at the CEP using HG and quasiparticle model. 



HG Models 



coordinates of CEP 



r tib tit 



n B /n T 



(T,Hb) 



(fm) (fm 3 ) (fm 3 ) 



Our HG Model 



(166, 155) 



0.8 0.104 0.54 



0.192 



Figure 8 demonstrates the variations of normalized energy density with respect to 
temperature at fJ,B = 400 MeV. We find a good agreement at higher temperatures in 
contradiction to the normalized pressure density case (see figure 7). Also one noticeable 
thing is that we found a discontinuity at the joining point of HG and QGP phase in variation 
of e/T 4 . However there is no sign of discontinuity present in lattice results. To make any 
clear statement one should wait for more refined lattice calculations. 

Figure 9 shows the variation of s/T 3 with respect to T at [Ib = 400 MeV. Lattice 
data points are also plotted to present a comparison with our model results. We find that 
at low as well as high temperatures our model results agree well with the lattice results. 
However, at intermediate temperature range i.e., from 160 MeV to 300 MeV, we find that 
our model results exceed over lattice results. We have also observed a discontinuity in our 
hybrid model results which clearly shows a first order phase transition at T = 160 MeV 
for jiB = 400 MeV in contradiction to lattice results which still shows a crossover type of 
behaviour at \xb = 400 MeV. 

Now we move towards showing the existence and nature of CP obtained in our hybrid 
model [13, 14]. We define a possible parameter x which gives the ratio of total baryon 
density rig to the total number density nj'(= rig + um) in the HG at the CP. Here Hm is 
the total meson number density. In table 1, we show the value of this ratio at CP obtained 
in our excluded volume prescriptions of HG. We find that the ratio x = 0.195 ±0.005 which 
signifies that at CP, the ratio um '■ ns is almost fixed as 4 : 1 . Beyond CP, the meson density 
increases and thus their fusion results into large cluster yielding a clear cross-over region. 
This demonstrates that the location of CP on the phase boundary is dependent on nature of 
the mesons and the baryons in our HG model since baryons possess a hard core volume but 
mesons can easily fuse into each other. Thus ub/um = 0.25 gives the location of CP. It still 
remains worth investigating why and how tib/tim = 0.25 defines the precise location of CP. 
In the excluded volume model, this ratio defines a critical fractional occupied volume Rq- 
As (Ib increases, the fractional occupied volume R by baryons increases and consequently 
the mobility of baryons in the hot and dense HG decreases fast and it finally results into a 
reduced pressure of HG so that the Gibbs' conditions of equilibrium phase transition become 
satisfied. The physical mechanism involved in this calculation is intuitively analogous to the 
percolation model where a first order phase transition results due to 'jamming' of baryons 
which thus restricts the mobility of baryons [23, 24]. However, in the percolation model 
we do not have any comparison to what we should get in the QGP picture. Here we use 
a similar picture and we explicitly and separately consider both the phases, i.e., HG as 
well as QGP and hence it gives a clear understanding how a first-order deconfining phase 
transition can be constructed in nature and finally we reach to an interesting finding that 
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the baryonic size is crucially responsible for the existence of CP on the phase boundary in 
such a construction. At low baryon density, overlapping mesons fuse into each other and 
form a large bag or cluster, whereas at high baryon density, hard-core repulsion among 
baryons, restricts the mobility of baryons. Consequently we consider two distinct limiting 
regimes of HG, one beyond CP is a meson-dominant regime and the other is a baryon 
dominant region. 

In conclusion, we have constructed a proper and realsitic hybrid model for the descrip- 
tion of strongly interacting QCD matter where we make use of an excluded volume model 
for HG and a thermodynamically consistent quasiparticle model for QGP phase. This 
model suitably works even in finite hb region and provides intuitive picture for crossover 
transition. We have also attempted to throw some light on the existence and the precise 
location of CP. We hope our results will be significant in future investigations regarding the 
EOS for QCD matter and the existence of CP on the QCD phase diagram. 
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